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A brief analysis of the dynamics of a Friedmann-Robertson- Walker uni- 
verse with a conformally coupled, real, self-interacting, massive scalar field, 
based on the Painleve theory of differential equations, is presented. Our re- 
sults complete earlier works done within the framework of Dynamical System 
Theory. 

We conclude that, in general, the system will not be integrable and that 
the chaos that has been found in a previous work, arises from the presence 
of movable logarithmic branch points in the solution in the complex plane of 
time. 
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Several works in the latter years have tried to prove the importance of non-linearities 
in General Relativity 0-0] that may lead to the onset of chaos in the early stages of the 
universe. The chaotic behavior is even observed to occur in quite simplistic models like 
vacuum Bianchi IX (see 0-0]), due to only the coupling of the different scale factors, as this 
model does not contain any matter fields. The focus has also been put on the inflationary 
phase of the early universe, and on the role that different fields might have played in this 
scenario PJ^J^. In many of these papers, the problem was attacked by means of numerical 
tools, based on Dynamical System Theory, which as we will show later in this brief report, do 
not explore some of the very interesting features of the system under study. Moreover, some 
people have shown (see for a discussion of the problem) that the detection of chaotic 
behavior is quite dependent on the choice of the coordinate system, and thus methods 
which provide gauge independent (or covariant) measures or results about the behavior of 
the system under study are very much needed while working within theories of Gravitation. 
Lately |jlO[ studies based on fractal basin boundaries have shown to be very useful, as no 



smooth coordinate transformation can undo a fractal pattern. On the other hand, from the 
analytical point of view, singularity studies can provide with very meaningful results, which 



*ahelmi@strw. LeidenUniv.nl 

^Member of the Carrera del Investigador de CONICET. 



1 



are independent of any particular gauge choice. And even more, they can be shown to be 



invariant under the group of homographic transformations of the dependent variables [|Tl 

In this report we study the dynamics of a conformally coupled field in a homogeneous 
and isotropic universe. We shall consider a massive, self-interacting, scalar field conformally 
coupled to a Friedmann- Robertson- Walker metric |T2[. The system is described by the 
following two degrees of freedom Hamiltonian: 

H = li-ipl + ka') + {pI + k^') + m'a^' + + ^a'] (1) 

where a is the scale factor, k is the curvature and A = ^ is related to the cosmo- 
logical constant. We have obviously considered a potential V^($) = |m^$^ + |A$^, where 
$ = (j)/{ay/v), X = X/v is the self-interacting coefficient {v is the conformal volume) and m 
is the mass of the field, can take both positive and negative values, what would take into 
account processes of spontaneous symmetry breaking resulting in phase transitions which 
may have occurred in the early stages of the universe [p!3| -[l6| . 



The Hamiltonian presented above (eq. (|I])) is well-known and has been thoroughly studied 
within lattice-dynamics, condensed matter theory, field theory, etc. (see references in [0). 
It is interesting to note that if we make the change of variable a = —ix, the Hamiltonian 



would correspond to that of a set of two quartically coupled non-linear oscillators [|T^], the 
study of which can be done analytically, giving a more complete view of the behavior of the 
system. 

The analytical method is based on the Painleve analysis of the differential equations, 
which states that if a solution is free from movable critical points other than poles, then it 
can be expressed as a generic Laurent series around one of its singularities. The test known 
as ARS algorithm , goes through the analysis of the leading order terms, the next order 



terms (or resonances) and the number of arbitrary constants that can be introduced in the 
expansion so as to make it generic. (It is important to note that the test gives only necessary 
conditions for the integrability of the system.) 

When one studies the dominant behaviors of the system derived from eq. (1) 

a + ka- m^acf)'^ - Aa^ = 0, (2) 
(f) + k(j) + m^a^(\) + X(\? = 0, (3) 

it is immediately possible to see that the leading order terms are independent of the value 
of the curvature k (see TABLE D, which is a very interesting and important result, as we 
shall interpret below. The cases for which the system passes P-test are very few. They 
correspond to special values of the parameters such as 



[ A = A and = -3A, 

1^ A = A and = —A, 

A = 16A and = — 6A, 

A = 8A and = -3A, 



Wk, (4) 
for k = 0. (5) 
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In reference the integrals of motion are given for these cases. 

For the rest of the cases, which are the majority , the system shows different types of 
movable critical points: algebraic branch points, with irrational or complex values for the 
exponents, or logarithmic branch points. (See TABLE Notice that, indeed, for a very few 
values of the parameters Branch 1 will have integer resonances. In general the resonances 
will be irrational or complex, thus leading to the nonintegrability of the system. The same 
would happen in Branch 2, if (r^ = 1 — 8m^/A < 0), or if is not the square of an integer 
or rational number. The same of course applies to Branch 3.) 

It is interesting to compare these results with those obtained by Blanco et al. |lT9| , pO| by 
means of numerical methods. They analyzed the same Hamiltonian, but restricted them- 
selves to the cases for which k (recall that /c = is a very important case indeed) 
and to a particular choice for the values of the parameters (2m^ + A + A = 0). The first 
restriction comes from a stability analysis which shows that the only fixed point for A; = is 
an elliptic point. This lead them to conclude, erroneously as we shall show, that a chaotic 
behavior cannot be reached in this fiat universe. The second condition is a consequence of 
the Hamiltonian constraint if = in General Relativity. With our approach, this constraint 
is always satisfied for each of the leading orders considered, without having to impose any 
restrictions on the values of the parameters. It is very easy to show with the method that we 
have been using, that the system passes the first two steps of the algorithm (leading order 
corresponds to Branch 1, the resonances are —1 and 4 with double multiplicity), but that the 
intented Laurent expansion will not have a sufficient number of arbitrary constants unless 
one introduces logarithmic terms in it. This fact is seen in the following: — Oq = b] 



|A-A| 



ta2 



G^sS +A2 -6A A) ( ~ '^0 ) , h = ai = 03 = 6i = 63 = 0, 04 is arbitrary, and 64 = 7(04) 



Thus a general solution would be of the form: x(t) = ia(t) = aoT"^+a2T+(a4+C4logr)r''+..., 
and a similar expansion for 0(t). 

Moreover, the chaos that has been observed by Blanco et al., can immediately be re- 
lated to the presence of these movable logarithmic branch points. It is well known pT| , p2 



that systems possessing this type of critical points display very interesting features in the 
complex plane of time. The singularities tend to accumulate in the complex t— plane, to 
form characteristic patterns such as chimneys (or spirals with increasing and decreasing radii 
for the case of complex resonances). It is also possible to show that the distance between 
neighboring singularities steadily decreases. Because of this, it is impossible to extend the 
solution further than a certain distance from a pair of initial singularities. Any path of an- 
alytic continuation appears to be trapped in a geometrically converging web of singularities 
that creates a natural boundary of the solution p3| , p^ . 

In synthesis, the above study has shown that there exists a relation between the oc- 
currence of movable branch points (i.e. multivaluedness, formation of self-similar patterns 
which eventually will be natural boundaries if they lie on the same Riemann sheet) and the 
onset of chaos in a system. We should emphasize again that this result is independent of 
the value of k, and includes the case of a fiat universe. This shows that there is no reason 
to believe that chaos will not be reached in this case. 
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Our results, based on the Painleve analysis, show that for most of the values the pa- 
rameters can take, the system is not integrable because of the presence of movable branch 
points, either algebraic (irrational and complex values for the resonances) or logarithmic. 
Besides, our results are independent of the curvature, as might have been expected, since it 
is well-known that the role of the curvature term can be neglected in Friedmann's equation 
during the inflationary phase. Thus, chaos is likely to be developed for all initial conditions 
of an homogeneous and isotropic universe, independently of the value of the curvature, cos- 
mological and coupling constants, except for a zero measure set of values. The cases for 
which a chaotic stage may have never been reached are only a very few p5[, that include 
for example when the coupling constant A is equal to the cosmological constant A. This 
fact would in some sense rule out the very well-known initial-value problem of the very 
early states of the universe, if the time scale for which the chaotic behavior is completed is 
comparable to the characteristic time for which the inflationary phase is expected to last. 
If during the final stages of A-dominated inflation the scalar field were to oscillate around 
the potential minimum in a chaotic way, the particles produced in this process could be 
responsible for the subsequent reheating of the universe p6|. This reheating in a chaotic 
fashion might be then, the necessary one to produce a sufficiently high temperature so as to 
restore the scenario for the standard baryogenesis to take place. 
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TABLE I. Leading order terms: x aor^, ^ ~ 6or^ (with t = t — to) 



p 


'7 "o 




Branch 1 —1 

Branch 2 — ^^-^ — 

Branch 3 —1 

lull ^ 1 1 -L 


—1 arbitrary 

l±Vl-8m2/A o/N 
2 /'^ 


^m^-AA 

— 2/A 

(All. k^xuxctiJ- y 


Branch 2 only appears if (m^ 
(m^ > and A > 0), or if (0 < - 


> and A > 0), or if (0 < -m^ < A). 
< A). 


And Branch 3, if 


TABLE IL 


Resonances: x{t) = aor^ + aiT^^^, etc. 




r 


Branch 1 x - | ± ^(m^-AA)^ 7m« 
Branch 2 
Branch 3 


- 16(A + A)m« - 18AAm4 + 16AA(A + A)m 
:fV1 - 8mVA, -1, 0, 4 
^Vl-8mVA, -1, 0, 4 


^ + 25A^A^ -1, 4. 
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